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A  study  was  made  of  the  short  period  detection  characteristics  of  LASA  within  the  codas  of  earthquakes 
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Kamchatka  would  have  gone  undetected  at  LASA  due  to  the  seismicity  of  the  latter  two  regions.  This  total 
mask  time,  under  various  assumptions,  is  quantified  into  separate  time  intervals  in  the  cases  where  the 
explosions  and  earthquakes  are  collocated.  A  result  of  this  quantification  is  that,  in  the  context  of  this 
experiment,  there  were  a  maximum  of  300  two  minute  time  intervals  during  the  decade  that  a  magnitude 
5.0  explosion  in  the  Kuriles  might  possihly  have  gone  undetected  at  LASA,  and  10  such  intervals  that  the 
explosion  most  probably  would  have  gone  undetected,  due  to  the  seismicity  of  the  Kurile  Islands  region. 
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ABSTRACT 


A  study  was  made  of  the  short  period  detection  characteristics  of  LASA  within 
the  codas  of  earthquakes  occurring  in  the  Kurile  Island  region  and  near  the  Kamchatka 
peninsula.  This  study  yielded  an  estimate  of  the  time  interval  that  LASA  would  not 
detect  an  explosion  in  central  Asia  following  an  earthquake  in  either  of  the  two  seismic 
regions;  this  estimate  is  given  as  a  function  of  earthquake -explosion  magnitude  difference. 
An  attempt  was  also  made  to  model  the  hypothetical  cases  where  the  explosions  were 
collocated  with  the  earthquakes  in  either  the  Kuriles  or  Kamchatka  regions  and  estimates 
of  the  mask  time  intervals  are  given  for  these  cases.  The  results  of  the  experiments 
are  extended  in  order  to  estimate  the  total  time,  during  the  decade  1961-1972,  that  an 
explosion  of  a  given  magnitude  in  either  central  Asia,  the  Kuriles,  or  Kamchatka 
would  have  gone  undetected  at  LASA  due  to  the  seismicity  of  the  latter  two  regions. 

This  total  mask  time,  under  various  assumptions,  is  quantified  into  separate  time 
intervals  in  the  cases  where  the  explosions  and  earthquakes  are  collocated.  A  result 
of  this  quantification  is  that,  in  the  context  of  this  experiment,  there  were  a  maximum 
of  300  two  minute  time  intervals  during  the  decade  that  a  magnitude  5.  0  explosion  in 
the  Kuriles  might  possibly  have  gone  undetected  at  LASA,  and  10  such  intervals  that 
the  explosion  most  probably  would  have  gone  undetected,  due  to  the  seismicity  of  the 
Kurile  Islands  region. 


Accepted  for  the  Air  Force 
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Acting  Chief,  Lincoln  Laboratory  Lia.son  Office 


iii 


INTRODUCTION 


Consideration  of  a  treaty  banning  the  testing  of  nuclear  weapons  underground 
raises  the  question  of  the  degree  of  deterioration,  du°  to  the  natural  seismicity  of  the 
earth,  of  the  detection  capability  of  a  seismic  station  or  network  of  stations.  Little,  if 
any,  formal  analysis  of  this  problem  exists  in  the  literature  of  seismology.  This 
report  describes  a  study  of  a  limited  aspect  of  this  topic.  It  attempts  to  answer  the 
question:  For  how  much  total  time  during  a  decade  would  an  explosion  of  a  given  size 
at  any  one  of  three  Asian  test  sites  go  undetected  at  the  Large  Aperture  Seismic  Array 
(LASA)  due  to  the  seismicity  of  the  Kurile  Islands  and  the  Kamchatka  Peninsula? 

No  attempt  was  made  here  to  solve  the  problem  of  explosion  detection  within 
earthquake  codas,  only  to  define  the  magnitude  of  the  problem.  Additionally,  no  attempt 
was  made  to  address  the  problem  of  identification  of  the  explosion  once  it  has  been 
detected.  An  effort  was  made  to  make  the  experiments  as  realistic  as  possible  in  that 
an  existing  detection  scheme  was  applied  to  actual  data. 

In  what  follows  the  general  approach  to  the  problem  is  described,  the  details 
o'  ^  method  are  given,  and  the  results  of  the  experiment  are  set  down.  A  section 
has  been  included  which  applies  the  results  of  this  work  to  the  problem  of  detecting 
clandestine  explosions  fired  just  after  earthquakes  for  the  purpose  of  detection  evasion. 
This  is  somewhat  an  awkward  topic  since  it  requires  discussion  of  detecting  an  act 
that  might  have  possibly  been  agreed  to  as  banned.  Nevertheless  the  quantification  of 
the  risks  of  undetected  violations  is  somewhat  relevant  to  the  discussion  of  a  nuclear 
test  ban  agreement  and  such  a  quantification  is  considered  here  in  a  very  limited  and 
qualified  sense.  The  chief  qualification  is  that  this  study  is  limited  to  th  use  of  one 
seismic  array,  LASA;  how  the  results  based  on  this  lone  station  will  scale  to  a  global 
network  is  not  considered  here. 


1 


Thus  the  methods  and  arguments  presented  below  may  merit  some  general  con¬ 
sideration,  the  numbers  that  come  out  of  these  processes  are  by  no  means  definitive 
in  the  sense  of  world  wide  seismicity  and  seismic  monitoring  networks. 
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METHOD 


In  this  section  the  method,  that  was  used  to  estimate  the  total  period  of  time 
that  an  explosion  of  a  given  magnitude  would  have  gone  undetected  at  LASA  due  to  the 
seismicity  of  a  given  region,  is  set  down.  Let  the  functions  M,  E,  and  t,  be  defined 
as  follows: 


M(m  ,T,s2) 


The  number  of  earthquakes  of  magnitude  m^  that 
occurred  in  the  region  s2  over  the  period  of  time 


E(mq-mx*  si»  s2} 


T(mx,T,  srs2) 


an  estimate  of  the  time  that  an  explosion  of  magni¬ 
tude  mx  at  site  s^  would  go  undetected  at  LASA  in 
the  P  wave  coda  of  an  earthquake  of  magnitude 
from  region  s2,  and 

an  estimate  of  the  total  mask  time  that  an  explosion 


of  magntiude  m  at  site  s,  would  have  gone  un- 

A  X 

detected  at  LASA  due  to  the  natural  seismicity  of 


region  s2  over  the  period  of  time  T. 

The  assertion  is  made  that  the  total  mask  time  t  may  be  represented  as 


T(mx»T»si,s2)  =  S  M(m  ,T,s2)E(m  -mx,s1,s2) 


wheie  the  summation  is  over  some  range  of  m  . 


Surely  at  some  low  value  of  mx 


an  explosion  of  tha:  magnitude  would  never  be  detected  due  to  the  ambient  seismic 


noise  in  the  earth.  All  explosions  at  some  higher  value  of  m  would  be  detected  if  it 
were  not  for  the  natural  seismicity  of  the  earth.  Equation  (1)  can  be  used 
to  estimate  the  effect  of  this  natural  seismicity  on  the  ability  of  LASA  to  detect  these 
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higher  magnitude  explosions. 


The  chief  difficulty  met  in  attempting  to  use  equation  (1)  lay  in  the  determination 
of  the  function  E.  This  function  is  called  a  failure  function  since  it  represents  an  esti¬ 
mate  of  the  amount  of  time  LASA  failed  to  detect  an  explosion,  in  the  coda  of  an  earth¬ 
quake,  given  an  earthquake -explosion  magnitude  deference.  Briefly,  the  scheme  used 
here  to  estimate  E  was  as  follows.  The  short  period  waveform  of  the  LASA  beam  from 
a  presumed  explosion  was  scaled  to  various  magnitudes  and  repeatedly  added  to  the  coda 
of  a  beam  formed  using  the  short  period  data  of  an  earthquake.  The  earthquake  beam 
was  always  aimed  at  the  presumed  explosion  site  in  order  to  enhance  the  explosion 
signal  with  respect  to  the  earthquake  coda.  Knowing  the  earthquake  magnitude  m^, 
the  scaled  explosion  magnitude  m^,  £tnd  the  number  of  scaled  explosions  repeatedly 
added  to  the  earthquake  coda  every  tQ  seconds;  the  number  (N)  of  these  explosions  that 
were  missed,  by  some  standard  detection  scheme,  was  determined.  By  repeating  the 
process,  using  other  earthquake  data  and  scaling  the  explosion  waveform  to  other 


magnitudes,  a  suite  of  measurements  of  N  were  made  as  a  function  of  m  -m  . 

q  x 


allowed  an  estimation  of  the  function  E  using 


This 


E(mq'mx'srs2)  =  to  *  N(mq'mx'srs2)  (2> 

Of  course  N(m^-mx,  s^,  S2)  was  not  found  to  be  constant  for  an  explosion  of 

magnitude  at  site  s^  and  all  earthquakes  of  magnitude  m  within  region  s9. 

Q  2 

Although  N  increased  with  the  difference  m  -mx,  the  scatter  about  any  specific 

value  of  the  difference  was  considerable.  This  was  probably  due  to  the  variance 

in  coda  shape  and  duration  between  earthquakes  from  the  same  region.  As  will 

be  seen  however,  the  maximum  number  of  explosions  missed,  N  ,  was  fairJv 

max  y 

well  de  ined  and  could  be  represented  by  an  expoential  function  as 


Nmax  =  6expf?(mq-mx)].  (3) 

Thus  an  estimate  of  the  maximum  detection  failure  time  at  LASA,  of  an  explosion  of 

magnitude  m  at  site  s.  due  to  an  earthquake  of  magnitude  m  at  site  s0,  is 
xi  q  z 

E  (m  -m  ,  s,,s„)  =  t  '  N  C4i 

max'  q  x  1  2'  o  max  '  > 

where  tQ  is  the  interval  between  the  explosion  waveforms.  In  practice,  the  equation  used 
to  estimate  the  maximum  total  time,  Tmax»  that  an  explosion  of  magnitude  mx  at  site 
Sj  would  have  gone  undetected  at  LASA  due  to  the  natural  seismicity  of  region  S2  over 
the  time  period  T  was 

Tmax(nVT>sl's2>  =  E  M(m.T,s2)t0  NJ«  (5) 


where  N  was  computed  from  (3). 

IIlclX 

In  this  study  an  attempt  was  made  to  estimate  t  (m  ,  T,  s , ,  s  J  for  the 

max  x  1  2 

following  four  cases. 


Case  I 


presumed  Soviet  test  site  (STS)  Semipalatinsk,  Kazakh 
Kurile  Islands  region 


Case  II 


Kamchatka  Peninsula  region 


Case  III 


hypothetical  test  site  in  the  Kurile  Islands 
Kurile  Islands  region 


Case  IV  Sj  -  Hypothetical  test  site  on  the  Kamchatka  Peninsula 

S2  -  Kamchatka  Peninsula  region. 

The  magnitude  range  of  m^  was  taken  to  be  4.  0  -  7.  1  and  the  time  T  to  be  the  decade 
1963-1972. 

Cases  I  and  II  are  not  contrived;  presumably  explosions  are  detonated  at 
the  Semipalatinsk  site  and  it  is  acceptable  to  add  these  signals  to  the  earthquake  codas 
in  order  to  effect  the  detection  experiments.  The  contrived  cases,  III  and  IV,  are 
justified  in  that  they  represent  the  extreme  instances  where  the  explosion  and  the 
earthquake  are  collocated.  Here  the  explosion  detection  beam  must  be  steered  directly 
at  the  interfering  earthquake  signal,  thus  minimum  earthquake  coda  attenuation  due  to 
beamforming  is  expected. 


DETAILS  OF  THE  EXPERIMENT 

The  following  event  was  taken  as  a  typical  explosion  from  STS;  the  parameters 
were  given  by  the  U.S.  Coast  and  Geodetic  Survey  (USCGS): 


Date 

5  September  1968 

Origin  time  - 

04h05m57. 4s 

Latitude 

49.  8°  North 

Longitude  - 

78. 1°  East 

Depth 

0. 0  kms  (fixed) 

Magnitude  - 

5.5  based  on  24  observations. 

The  appropriate  LASA  beam  trace  was  formed  and  this  trace  was  scaled  to  ten  other 
magnitudes  (mx)  through  division  by  the  following  factors  (f): 


f 

mx 

f 

mx 

1 

5.5 

20. 

4.2 

2 

5.2 

32. 

4.0 

3.2 

5.0 

50. 

3.8 

5.0 

4.  8 

80. 

3.6 

8.0 

4.6 

125. 

3.4 

12.5 

4.4 

Before  scaling  the  original  explosion  beam  trace  was  filtered  using  a  third  order  band 
pass  filter  with  3  db  corners  at  .  9  and  1.4  Hz.  The  filter  is  used  routinely  in  the  LASA 
detection  scheme.  .  Figure  1  shows  the  unfiltered  and  filtered  versions  of  the  explosion 
beam  trace,  sampled  10  times  a  second  and  repeated  every  20  seconds. 

The  array  data  for  the  various  earthquakes  was  then  phased  using  delays 
appropriate  for  the  proposed  explosion  site  and  summed.  This  earthquake  beam  was 
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-2-11008 


Fig.  1.  The  traces  of  the  short  period  LASA  beam  showing  the  P  wave  from  a 
presumed  explosion  in  eastern  Kazakh,  S.  S.R.  The  trace  is  repeated  every 
20  seconds  as  are  the  vertical  dotted  lines.  The  lover  trace  resulted  in  band¬ 
pass  filtering  the  unfiltered  upper  trace.  Scaled  versions  of  the  lower  trace 
were  added  to  filtered  earthquake  codas  in  order  to  perform  the  detection  ex¬ 
periments  described  in  the  text. 
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filtered  using  the  same  filter  describe  ',  in  the  last  paragraph.  The  filtered  earthquake 
beam  was  then  added  to  the  repeated  filtered  explosion  beam  trace  for  use  in  the  detec¬ 
tion  study. 

In  Cases  I  and  II  this  procedure  is  straightforward  and  somewhat  realistic. 
Cases  III  and  IV  are  less  so  ;n  that  there  are  no  known  large  underground  explosions 
in  the  Kuriles  or  Kamchatka  region.  In  these  cases  the  same  STS  explosion  beams  were 
used  ■  3  in  Cases  I  and  II,  however  the  earthquake  beams  were  formed  on  hypothetical 
sites  in  the  Kuriles  and  Kamchatka  regions.  This  is  equivalent  to  asserting  that,  at 
LASA,  the  short  period  P  waves  from  an  explosion  in  far,  northeastern  Asia  will  be 
similar  to  one  in  central  Asia.  Because  of  geologic  contrasts  between  the  actual  and 
hypothetical  explosion  sites  this  assertion  is  questionable  but  the  ruse  is  necessary  to 
form  such  experiments.  The  locations  of  the  explosions  and  earthquakes  for  the  various 
cases  are  summarized  below: 


Explosion 

Earthquakes 

Case  I 

50° N  78° E 

40-50°N,  145-155°E 

Case  II 

50° N  78° E 

50-60°N,  155-170°E 

Case  III 

45° N  150° E 

40-50°N,  145-155°E 

Case  IV 

55°  N  160°E 

50-60°N,  155-170°E 

The  explosions  were  all  taken  to  be  at  the  surface  and  no  subdivision  of  the  earthquake 

population  with  respect  to  depth  was  attempted.  A  complete  list  of  the  earthquakes 

used  to  estimate  the  function  N(m  -mx,  sJf  s2)  is  given  in  the  Appendix.  The  earthquake 

regions  given  above  are  also  those  used  to  compute  the  seismicity  function  M(m  ,T,  sj 

q  Z 

from  USCGS  seismicity  lists  for  the  decade  1963-1972. 
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The  detection  scheme  now  remains  to  be  described.  Tile  scheme  was  modeled 
after  that  used  at  tile  Seismic  Array  Analysis  Center  (SAAC)  and  is  described  in  the 
reference.  Basically  the  method  compares  the  ratio  of  a  short  term  rectified  signal 
average  to  a  similar  long  term  average  and  declares  a  detection  when  then  ratio  exceeds 
a  certain  threshold  for  a  given  number  of  comparisons.  No  search  was  made  to  find 

an  optimum  coda  detection  process;  a  standard  process  was  used  with  threshold  para- 
meters  that  gave  relatively  few  false  alarms. 

The  short  term  average  (STA)  was  computed  using  the  following  formula  from 
reference  1, 

1  P'1 

f%T  A  A  t-  \  —  T1  I  nr  /„  v  »  —  I 

(6) 


STA(nAt)  =  p  Z  |B[(n-p)At]|. 

p=0 


Here  B[nAt]  represents  the  filtered  beam  sampled  at  nAt.  In  practice  At  was  set  at 

.  2  see  and  P=3,  i.  e.  the  STA  was  computed  every  .  6  sec.  The  long  term  average  (LTA) 

was  computed  using  the  current  STA  value,  STAn,  and  the  last  LTA  value,  LTA^,  and 
the  recursive  relation 


LTAn  =  2^STAn  +  (l-2-°)LTAn_r  (7) 

In  practice  values  for  n  and  a  were  1.  and  4.  respectively. 

Because  the  effective  length  of  this  LTA  was  greater  than  the  20  sec  interval 
at  which  the  explosions  were  added,  the  LTA  was  computed  using  the  filtered  earthquake 
beam  to  which  no  explosion  traces  had  been  added.  This  means  that  the  added  explosions 
did  not  contaminate  the  LTA  upon  which  the  detections  were  based.  Eleven  STA 
values  were  computed  for  the  respective  scaled  explosion -earthquake  sum  traces  and 
also,  to  check  false  alarms,  on  the  earthquake  trace  alone.  Thus,  12  STA/LTA  ratios 


were  computed  every  .6  sec,  one  on  the  earthquake  trace  and  one  on  each  of  the  11 
earthquake  traces  to  which  explosions  scaled  to  11  different  magnitudes  were  added  every 
20  seconds.  A  detection  was  declared  on  any  of  the  twelve  traces  when  the  STA/LTA 
ratio  exceeded  10  db  for  1.8  seconds,  i.e.  for  three  consecutive  STA/LTA  measurements. 

An  example  of  the  method  and  the  performance  of  the  detection  scheme  is 
shown  by  Figure  2  where  earthquake  number  15  is  shown  under  the  conditions  of  Case 
III.  The  top  three  traces  are  the  earthquake  trace  plotted  at  different  gain  levels. 

Traces  4  and  5,  with  explosions  equivalent  to  magnitudes  5.  5  and  5.  2  added,  are  plotted 
at  the  same  gain  as  trace  1.  Traces  6  and  7,  explosion  magnitudes  5.  0  and  4.  8,  are 
plotted  at  the  gain  of  trace  2  and  traces  8  and  9,  explosion  magnitudes  4.  6  and  4.4,  at 
the  gain  of  trace  3.  Explosion  declared  detections  are  marked  with  an  "X",  undetected 
explosions  are  marked  with  an  "O".  This  writer  feels  that,  in  the  light  of  such  examples, 
the  detection  criteria  used  were  reasonable  ones  however  the  reader  must  make  his  own 
judgement.  One  other  point  bears  mention  at  this  stage.  As  it  turned  out  in  most  runs 
of  the  experiment  the  explosions  were  added  such  that  one  usually  fell  a  few  seconds 
before  the  beginning  of  the  earthquake  arrival.  This  will  tend  to  minimize  the  number 
of  explosions  missed  but  this  feature  should  be  balanced  when  the  failure  function  is 
selected  to  represent  the  maximum  number  of  explosions  undetected. 

The  remaining  details  of  the  experiment  are  simply  described.  The  detection 
criterion  was  applied  to  the  earthquake  and  the  11  earthquake  plus  explosion  traces.  The 
number  of  false  alarms  was  recorded  from  the  former,  the  number  of  explosions  missed 
from  the  latter.  If  Tie  number  of  detected  explosions  fell  below  the  number  of  explosions 
added  before  the  earthquake  signal,  due  to  high  noise  conditions,  that  trace  was  dis¬ 
regarded.  Otherwise  the  number  of  explosions  missed  was  plotted  against  the  magnitude 
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difference  between  the  earthquake  and  the  scaled  explosion  and  the  failure  function  for 
each  case  was  established. 
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RESULTS  OF  THE  EXPERIMENT 


The  number  of  explosions  undetected  by  the  method  just  described  are  plotted 
versus  the  earthquake -explosion  difference  in  Figures  3-6  for  Cases  I-IV  respectively. 
The  straight  lines  drawn  on  these  figures  represent  the  author’s  estimate  of  the  maximum 
lumber  of  explosions  that  would  be  missed  for  a  given  magnitude  difference  in  each  case. 
The  parameters  of  these  lines,  which  are  in  the  form  of  equation  (3),  are 


Case 

6 

§ 

I 

1.9 

1.66 

II 

2.9 

1.64 

III 

4.5 

1.20 

IV 

5.5 

1.26. 

The  distribution  of  points  and  the  parameters  of  the  straight  lines  on  these 
figures  merit  some  discussion.  Although  Cases  I  and  II  (Figures  3  and  4)  are  similar 
in  that  the  explosions  are  from  the  same  site  and  the  earthquakes  are  from  regions 
contiguous  to  each  other,  the  scatter  of  the  points  is  greater  in  Case  II  than  in  I.  The 
earthquakes  used  in  Case  I  (Kurile  Islands)  all  occurred  within  a  two  month  period  at 
nearly  the  same  location  (43N  147E)  while  those  of  Case  II  (Kamchatka)  are  more  widely 
dispersed  in  time  and  space.  Thus,  although  the  lines  on  the  two  figures  have  about  the 
same  slope  (?)  the  zero  intercepts  (6  )  vary  by  1.  0.  The  same  observation  arises  when 
the  maximum  failure  lines  on  Figures  5  and  6  are  compared.  It  is  also  of  interest  to 
compare  the  trend  of  the  points  in  Cases  I  and  II  with  that  of  III  and  IV.  The  former 
appear  to  exhibit  a  steeper  slope  which  is  reflected  in  a  higher  value  of  ?  in  the  maximum 
failure  lines  drawn  for  Cases  I  and  II  than  in  III  and  IV.  The  slope  or  trend  of  these 
lines  is  largely  controlled  by  the  number  of  missed  detections  at  high  magnitude  dif- 
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Fig.  3.  The  points  represent  the  number  of  STS  explosions  undetected 
on  the  LASA  beam,  aimed  at  STS,  due  to  Kurile  Islands  earthquakes. 
These  points  are  plotted  versus  the  magnitude  difference  (mq  —  mx) 
between  the  earthquake  (mq)  and  the  explosion  (mx).  This  figure  is 
relevant  to  Case  I  described  in  the  text.  The  straight  line  represents 
the  author's  estimate  of  the  maximum  number  of  explosions  missed 
as  a  function  of  magnitude  difference.  The  explosions  were  added  at 
20  second  intervals. 


ferences,  determined  far  back  in  the  coda  of  the  larger  earthquakes.  This  leads  to  the 
conclusion  that  the  coda  from  large  events  becomes  increasingly  incoherent  with  in¬ 
creasing  time  after  the  first  arrival.  Alternatively  stated,  the  effectiveness  of  beam- 
forming  on  location  A  to  suppress  the  coda  of  an  event  from  location  B  decreases  with 
time  after  the  arrival  of  the  initial  P  wave  from  location  B. 

An  attempt  was  now  made  to  extend  these  results  to  estimate  the  maximum 
total  time  that  an  explosion  of  a  given  magnitude  would  have  gone  undetected  at  LASA 
due  to  the  long  term  seismicity  of  two  regions  considered.  This  was  done  through  equa¬ 
tion  (5).  Before  discussing  the  results  of  these  calculations  however,  several  qualify¬ 
ing  comments  should  be  made  concerning  their  interpretation.  The  use  of  equation  (5) 
requires  an  accurate  knowledge  of  the  seismicity  of  a  region.  Here  USCGS  lists  have 
been  used  simply  because  they  appear  to  be  themost  consistent  source  of  seismicity 
information.  Figures  7  and  8  show  a  cummulative  number  of  earthquakes  as  a  function 
of  magnitude  over  the  decade  1963-1972  for  the  Kurile  Islands  and  Kamchatka  regions 
respectively.  It  is  clear  from  the  flattening  of  these  seismicity  curves  near  magnitude 
4. 0  that  this  magnitude  approximates  the  lower  limit  of  sensitivity  for  the  network  of 
stations  used  by  the  USCGS  to  locate  and  assign  magnitudes  to  seismic  events  from  these 
regions.  At  what  magnitude  any  flattening  due  to  insensitivity  and  not  nature  begins  , 
is  the  crucial  yet  undeterminable  point.  In  Figure  7  a  slight  bend  in  the  cummulative 
curve  is  seen  at  about  magnitude  5.4,  this  seems  too  high  to  be  due  to  network  insensi- 
ti  /ity  which  might  be  expected  to  become  significant  below  magnitude  5. 0.  In  any  case, 

this  effect  must  qualify  any  conclusion  based  on  equation  (1)  at  lower  explosion  magni¬ 
tudes. 

There  exists  a  more  specific  qualification  that  may  be  placed  on  any  result 


Fig.  7,  The  points  represent  the  total  number  of  earthquakes 
of  a  given  magnitude  (mq)  or  greater  reported  by  the  USCGS 
during  1963-1972  as  occurring  in  the  Kurile  Island  region 
40-  50°N,  145-  1 55°E. 


Fig.  8.  The  points  represent  the  total  number  of  earthquakes 
of  a  given  magnitude  (rriq)  or  greater  reported  by  the  USCGS 
during  the  1963-  1972  as  occurring  in  the  Kamchatka  region 
50- 60°N,  155- 170°E. 
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gained  from  the  use  of  equation  (5).  Because  of  the  form  given  the  failure  function  by 
equation  (3),  the  values  assigned  to  its  parameters,  and  the  detection  criteria  used;  an 
explosion  may  be  hidden  by  an  eartnq’’ake  of  smaller  magnitude.  Thus,  in  order  to 
estimate  the  total  time  that  a  magnitude  4.  9  explosion  would  be  hidden,  the  total  number 
of  magnitude  3.8  earthquakes  is  required.  S  nee  earthquakes  of  magnitude  3.8  are  not 
considered,  any  estimate  of  t  in  this  case  mu.  t  be  low.  The  severity  of  this  effect 
increases  with  the  value  of  the  parameter  6.  In  the  summation  of  equation  (5)  all 
fractional  detections  have  been  rounded  down  to  the  nearest  integer.  This  is  a  fairly 
safe  procedure  because  all  of  the  points  of  Figures  3-6  lie  at  least  one  unu  oelow  the 
maximum  failure  lines.  In  the  light  of  the  above  discussion  the  condition  that 

6  exp[5(mx-m^)]  <  1.0  (8) 

when  estimating  t  max ,  where  nT  is  the  magnitude  above  which  M(mq,  T,  s2)  is  believed 
to  be  complete,  assures  the  validity  of  the  estimate  t  in  the  context  of  equation  (5). 

Figures  9-12  show  the  result*  of  applying  the  failure  function  for  Cases  I-IV 
to  the  appropriate  10  year  seismicity  data.  The  ordinate  in  these  figures  is  the  estimate 
of  the  maximum  total  time  in  minutes  during  the  10  year  period  that  an  explosion  of 
magnitude  m^,  at  the  site  indicated,  would  have  gone  undetected  due  to  the  seismicity 
of  either  the  Kurile  Islands  or  the  Kamchatka  Peninsula.  The  arrows  on  these  figures 
indicate  an  mx  below  which  the  estimates  are  not  strictly  valid  because  of  the  problems 
discussed  in  the  last  two  paragraphs.  The  placing  of  this  arrow  is  based  on  the  condi - 
..on  of  equation  (8)  where  nT  =  4.5.  This  implies  that  the  seismicity  data  of  Figures  7 
and  8  is  considered  complete  down  to  mq  =  4.  5.  If  the  reader  chooses  to  pick  an  alter¬ 
nate  value  for  mq  the  arrows  on  Figures  9-12  will  simply  be  shifted  to  the  right  or  left 
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Fig.  9.  The  result  of  applying  Eq.  (1)  to  the  data  of  Figs.  3 
and  7  (Case  I).  The  points  represent  an  estimate  of  the 
time  (in  minutes)  that  a  magnitude  Mx  explosion  at  STS 
would  have  gone  undetected  by  LASA  due  to  Kurile  Island 
seismicity  during  the  decade  1963-1972.  The  arrow  rep¬ 
resents  the  magnitude  below  which  the  results  are  not  con¬ 
sidered  strictly  valid.  The  line  represents  the  approxi¬ 
mate  linear  extension  of  the  points  at  magnitudes  greater 
than  that  of  the  arrow. 
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Fig.  10.  A  diagram  analogous  to  Fig.  9  and  relevant  to  Case  II. 
The  points  represent  an  estimate  of  the  time  (in  minutes)  that 
a  magnitude  mx  explosion  would  have  gone  undetected  by  LASA 
due  to  Kamchatka  seismicity  during  the  decade  1 96 3-1 972 . 
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Fig.  11.  A  diagram  analogous  to  Fig.  9  and  relevant  to  Case  III. 
The  points  represent  an  estimate  of  the  time  (in  minutes)  that 
a  hypothetical,  magnitude  mx  explosion  in  the  Kurile  Islands 
would  have  gone  undetected  by*  LASA  due  to  Kurile  Islands  seis¬ 
micity  during  the  decade  1963-  1972. 
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Fig.  12.  A  diagram  analogous  to  Fig.  9  and  relevant  to  Case  IV. 
The  points  represent  an  estimate  of  the  time  (in  minutes)  that 
a  hypothetical  explosion  of  magnitude  mx  near  the  Kamchatka 
Peninsula  would  have  gone  undetected  byLASA  due  to  Kamchatka 
seismicity  during  the  decade  1963-  1972. 


the  number  of  magnitude  units  that  the  new  nF  differs  from  4. 5.  The  straight  solid 
lines  on  Figures  9-12  represent  the  author’s  approximation  of  the  linear  extention  of 
the  valid  points  to  the  right  of  the  arrows  to  mx  values  to  the  left  of  the  arrows  where 
the  points  shown  are  surely  low. 

As  an  example,  these  figures  allow  us  to  estimate  the  maximum  total  time 
that  an  explosion  of  say  magnitude  5. 5  would  have  goae  undetected  at  LASA  during  the 
ten  year  period  for  each  of  the  four  cases.  These  results  are  summarized  below: 

Case  Total  mask  time 

I  205  min  (3.4  hrs) 

II  135  min  (2.  3  hrs) 

HI  approx  2000  min  (33.4  hrs) 

IV  approx  100  min  (16.  7  hrs) 

Again  these  are  estimates  of  themaximum  total  mask  time  at  LASA  relevant  only  to 
the  cases  indicated.  The  maximum  total  mask  time  due  to  world  wide  seismicity  would 
of  course  be  greater. 

The  application  of  the  detection  scheme  to  the  earthquake  trace  alone  gave  rise 
to  nine  false  alarms  in  all  of  earthquake  codas  processed.  Of  the  spurious  detections, 
five  were  due  to  the  surface  reflection  phase  pP,  two  were  due  to  a  small  interfering 
event,  and  the  remaining  two  were  genuine  false  alarms.  This  means  that  twice  the 
nature  of  the  earthquake  coda  was  such  that  the  detection  method  mistook  it  for  a 
discrete,  second  seismic  event. 


EVASION 


The  results  shown  above  permit  some  discussion  of  the  problem  of  detecting 
an  explosion  deliberately  hidden  in  an  earthquake  coda  for  the  purpose  of  evasion.  It 
is  clear  fiom  Figures  9-12  that,  in  the  case  of  LASA,  the  chances  of  success  of  such 
a  scheme  will  increase  as  the  distance  between  the  explosion  site  and  the  masking  earth¬ 
quake  location  decreases.  This  is  due  to  the  earthquake  coda  attenuation  on  the  array 
beam  when  steered  at  an  explosion  in  a  different  region.  In  the  cases  studied  here,  it 
was  found  that  when  the  explosion  and  earthquake  locations,  at  about  the  same  distance, 
weie  separated  in  azimuth  by  about  50°  with  respect  to  the  array;  the  total  explosion 
mask  time  was  about  an  order  of  magnitude  less  than  when  the  two  source  types  were 
collocated.  If  one  considers  a  global  network  of  stations  the  chances  of  successfully 
hiding  an  explosion  in  the  coda  of  a  distant  earthquake  become  even  less.  In  order  to 
avoid  detection  on  a  global  network  the  clandestine  test  would  have  to  be  fired  at  a  time 
such  that  the  explosion  P  wave  would  not  arrive  before  the  earthquake  P  wave  at  any  point 
on  the  surface  of  the  earth.  The  way  to  ensure  this  would  be  to  detonate  the  explosion 
upon  arrival  of  the  earthquake  P  wave  at  the  explosion  test  site.  If  the  masking  earth¬ 
quake  were  at  say  50  distance  from  the  explosion  site,  this  would  require  waiting  some 
9  minutes  after  the  origin  time  of  die  earthquake.  At  an  array  at  teleseismic  and  equal 
distances  from  both  explosion  and  earthquake  sites  the  explosion  P  wave  would  arrive 
9  minutes  into  the  earthquake  coda.  Obviously  this  effect  greatly  reduces  the  chances 
of  success  of  such  a  scheme. 

Thus  any  potential  evasion  scheme  would  be  more  likely  to  succeed  if  the 
clandestine  explosion  was  detonated  within  a  region  of  high  seismicity.  Such  a  situation 
is  considered  here  in  Cases  III  and  IV.  In  order  to  hide  an  explosion  in  a  local  earth- 
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quake  one  would  require  the  means  to  detect,  locate,  and  fire  the  explosion  within  a  given 
time  interval.  Of  course,  the  larger  the  earthquake  the  larger  the  explosion  that  could 
be  hidden  or,  for  a  given  explosion  size,  the  longer  the  time  interval  available  for  seismic 
and  engineering  pre -firing  routines.  Thus,  it  is  important  in  any  evasion  study  not  only 
to  know  the  total  time  that  an  explosion  of  a  given  magnitude  might  be  massed  but  also 
the  quantization  of  this  time.  The  proper  quantization  is  considered  to  be  the  number 
of  time  intervals  of  a  given  length  available,  over  an  extended  period  of  time,  for  the 
hiding  of  an  explosion  of  a  given  magnitude.  An  attempt  toward  this  quantization,  using 
LASA  as  the  sole  detector,  has  been  made  for  Cases  III  and  IV  based  on  the  detection 
data  of  Figures  5  and  6  and  the  seismicity  data  of  Figures  7  and  8. 

There  is  one  point  that  must  be  considered  at  this  time  before  such  an  attempt 
is  made.  The  upper  lines  on  Figures  5  ar  d  6  represent,  at  least,  the  maximum  number 
of  explosions  missed  in  the  detection  experiments.  In  every  case  fewer  explosions  went 
undetected  at  LASA  than  these  lines  indicate.  An  evasion  scheme,  using  firing  time 
intervals  based  on  these  lines,  would  have  a  high  risk  of  detection  at  LASA  since  LASA 
always  did  better  than  the  estimates  these  lines  represent.  To  ensure  success,  with 
respect  to  LASA,  a  potential  evasion  scheme  must  be  based  on  the  minimum  number  of 
explosions  missed  at  LASA  for  a  given  earthquake-explosion  magnitude  difference.  The 
question  is  not  then:  What  were  the  number  of  t  minute  time  intervals  that  an  explosion 
of  magnitude  mx  in  the  Kurile  Islands  might  have  gone  undetected  by  LASA?;  but  rather: 
What  are  the  number  of  these  intervals  that  this  explosion  surely  would  have  gone  un¬ 
detected  at  LASA.  In  order  to  estimate  the  number  arid  length  of  these  latter  intervals 
using  Figures  5  and  6,  the  lower  lines  that  represent  the  minimum  number  of  explosions 
missed  must  be  used. 
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These  minimum  failure  lines  have  been  drawn  parallel  to  the  maximum 
failure  lines  but  shifted  down  by  about  a  factor  of  4.  The  parameters  of  these  minimum 
failure  lines  are 


Case 

6 

P 

III 

1.  1 

1.20 

IV 

1.  1 

1.26 

i  lie  pairs  of  these  lines  have  been  used  to  estimate  the  number  of  separate  time  inter¬ 
vals  of  a  given  length  (taken  to  the  nearest  minute)  during  1963-1972  than  a  magnitude 
5.  0  explosion  near  either  Kamchatka  or  the  Kurile  Islands  might  possibly  (based  on  the 
maximum  failure  line)  and  most  probably  would  (based  on  the  minimum  failure  line)  have 
gone  undetected  at  LAS  A  due  to  the  seismicity  of  the  two  regions.  The  results  are 
shown  m  Figures  13  and  14.  Due  to  Kurile  Islands  earthquakes  during  the  decade  there 
were  approximately  300  two  minute  intervals  within  which  a  magnitude  5.  0  explosion  in 
the  Kuriles  might  have  gone  undetected  at  LASA,  however  if  the  explosion  had  been  fired 
within  etch  interval  many  would  have  been  detected.  There  were  approximately  10  two 
minute  intervals  within  which  the  same  explosion  would  surely  have  gone  undetected  at 

LASA.  The  respective  number  of  two  minute  intervals  for  the  Kamchatka  region  are 
approximately  130  and  5. 

From  the  evader’s  point  of  view  there  wculd  have  existed  approximately  one 
two  minute  interval  a  year  during  which  a  magnitude  5.  0  clandestine  explosion  could  have 
been  fired  in  the  Kuriles  with  minimum  risk  of  detection  by  LASA.  From  the  detector’s 
point  of  view,  there  would  have  been  approximately  30  two  minute  intervals  a  year,  in 
the  codas  of  Kurile  earthquakes,  which  would  have  to  have  been  searched  for  a  clandestine 
explosion  of  magnitude  5.  0.  Of  these  latter  intervals  there  would  have  been  approximately 
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Fig.  13.  Estimates  of  the  number  of  separate  time  intervals  of  a  given 
length  that  a  hypothetical  magnitude  5.0  explosion  in  the  Kurile  Islands 
would  have  gone  undetected  at  LASA  due  to  the  seismicity  of  the  Kurile 
Islands  during  the  decade  1963-1972.  The  hollow  bars  represent  the 
maximum  number  of  intervals  during  which  the  explosion  might  pos¬ 
sibly  have  gone  undetected,  based  on  the  maximum  number  of  missed 
detections  at  LASA.  The  shaded  bars  represent  the  number  of  inter¬ 
vals  that  the  explosion  most  probably  would  have  gone  undetected  at 
LASA,  based  on  a  minimum  number  of  missed  detections  at  LASA. 
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Fig.  14.  A  diagram  analogous  to  Fig.  14  showing  the  estimates 
of  the  number  of  separate  time  intervals  of  a  given  length  that 
a  hypothetical  magnitude  5.0  explosion -On  the  Kamchatka  Pen¬ 
insula  would  have  gone  undetected  at  LASA  due  to  the  seismicity 
of  the  Kamchatka  region  during  the  decade  1963-1972. 


one  where  no  positive  determination,  of  whether  an  explosion  had  or  had  not  occurred, 
could  be  made. 

No  estimate  can  be  made  of  the  maximum  or  possible  number  of  one  minute 
intervals  because  of  the  condition  expressed  by  equation  (5).  The  number  of  one  minute 
intervals  within  which  it  would  have  been  impossible  to  determine  whether  a  5.0  test 
had  occurred  or  not  is  approximately  300  for  the  Kurile  Islands  and  80  for  the  Kamchatka 
region.  It  is  pointed  out  that  the  P  wave  travel  time  across  these  regions  is  approximately 
two  minutes. 
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CONCLUSION 


A  method  has  been  developed  which  yields  an  extimate  of  the  maximum  total 
time,  during  a  specific  period  of  time  at  a  seismic  station,  that  earthquakes  in  a  given 
region  will  mask  the  detection  of  a  second  event,  here  taken  to  be  an  explosion,  from 
another  region.  The  maximum  mask  time  is  given  as  a  function  of  explosion  magnitude 
in  each  of  four  cases.  In  two  of  these  cases  the  explosion  site  was  taken  to  be  in  central 
Asia  while  the  interferring  earthquakes  came  from  the  Kurile  Islands  and  the  Kamchatka 
Peninsula,  hi  the  two  other  cases  the  earthquakes  and  explosions  were  collocated  in  the 
Kuriles  and  near  Kamchatka.  The  total  mask  time  at  LASA  in  the  first  two  cases  was 
found  to  be  about  an  order  of  magnitude  less  than  the  latter  two.  In  these  latter  cases, 
it  was  found  that  a  magnitude  5.  5  explosion  in  the  Kuriles  would  have  gone  undetected 
for  a  maximum  total  time  of  about  33  hours  due  to  Kurile  Islands  seismicity  during  the 
decade  1963-1972.  The  same  explosion  on  Kamchatka  would  have  been  masked  for  a 
maximum  total  of  about  17  hours  during  the  same  period  due  to  the  Kamchatka  region 
seismicity. 

The  results  of  these  experiments  are  discussed  in  tne  light  of  a  scheme  which 
would  attempt  to  evade  detection  of  a  clandestine  explosion  by  firing  just  after  an  earth¬ 
quake.  In  such  a  discussion  it  is  important  to  not  only  know  the  total  amount  of  time 
available  for  such  a  ruse  but  also  how  this  time  is  quantified  and  what  the  risks  of  de¬ 
tection  are  for  certain  earthquake -explosion  magnitude  differences.  It  was  found  that  due 
to  Kurile  Islands  earthquakes  during  the  decade  there  were  a  maximum  of  approximately 
300  two  minute  intervals  within  which  a  magnitude  5.  0  explos’on  in  the  Kuriles  might 
possibly  have  gone  undetected  at  LASA.  During  approximately  ten  of  these  intervals 
the  explosion  surely  would  have  gone  undetected  at  LASA  due  to  Kurile  Islands  seismicity. 
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The  respective  number  of  two  minute  intervals  for  the  same  explosion  collocated  with 
earthquakes  in  the  Kamchatka  region  are  approximately  130  and  5. 

The  numbers  presented  in  this  paper  should  not  be  considered  as  definitive 
since  they  are  based  on  a  rather  crude  detection  scheme  applied  at  only  one  array. 
More  sophisticated  detection  schemes  applied  simultaneously  at  a  network  of  stations 
or  arrays  should  lower  these  estimates.  The  method  of  estimating  the  total  mask  time 
however  may  have  general  applicability. 
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